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1. Introduction 

Two-dimensional gravity has undergone a thorough examination over the last few 
years, especially with the emergence of efficient calculational techniques stemming from 
the matrix model approach. Two dimensions is an arena where rather complex phenom- 
ena (confinement, chiral symmetry breaking, integrability, nonperturbaive phenomena, 
solitons; the list is extensive) can be stripped of complications of higher dimensional kine- 
matics and dynamics while hopefully retaining many of the physical features of the problem 
of interest. Thus we study two dimensional gravity as a model for exploring the structure 
and formalism of quantum gravity: the wavefunction of the universe, the Wheeler-DeWitt 
equation the meaning of measurement in quantum gravity, the statistical properties of 
the metric and matter in a fluctuating geometry, etc. Indeed the collection of solvable 2d 
gravity-matter systems provides a rich laboratory for the investigation of these issues. In 
addition, since unified strings are by definition coordinate invariant 2d quantum field the- 
ories these systems are equally well regarded as solvable models of string theory, where one 
might begin to formulate a useful string field theoryf^, study nonperturbative effects]^, 
strong coupling phenomena, etc. 

There are several ways one might approach the continuum theory of 2d gravity. One 
is to write down a field theory on 2d metrics and matter [Q, regulate it covariantly, and 
try to find a consistent renormalization to the continuum limit [|5|0. Another is to dis- 
cretize the theory and study the fiuctuations of the discrete geometry 0, then try to take 
the continuum limit of the random lattice model (a method that has also been used for 
strings embedded in higher dimensional spacetimesQ). Both should yield the same results 
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if indeed 2d gravity is universal, i.e. if there is a second order phase transition for the fluc- 
tuating surfaces in some region of the space of coupling constants (cosmological constant, 
topological coupling, boundary cosmological constant, etc.). The former line of attack ad- 
heres more closely to the standard conceptual framework of gravity, yet in two dimensions 
the latter formulation has yielded more success to date. This success comes about largely 
because in two dimensions it is not difficult to enumerate the lattice configurations, which 
are gauge invariant. 

These lecture notes introduce the present understanding of 2d gravity /solvable string 
models, beginning in section 2 with a review of the continuum formulation in terms of 
the functional integral over metrics. While there is little that can be computed exactly, 
we argue that much of the classical structure should persist upon quantization; and that 
the familiar structures of rational conformal field theory might appear, albeit in a much 
more subtle and complicated way. In section 3 we discuss the behavior of the theory 
under rescalings of the 2d metric]^ as well as the related dynamics in the zero mode 
sector |TT[] |]T^. In section 4 we switch to the discrete formulation. Here the partition 
function of random surfaces can be recast as and integral over NxN matrices where in 
the continuum limit oo[jl^. Both the lattice[0 and continuum ||13|| [1T5[| ||16|| theories 



exhibit a rich integrability structure. The scaling behavior predicted by the continuum 



approach is recovered as well as the zero mode dynamics |T0] [12 



2. The effective action for gravity and its quantisation 

A crucial feature of 2d quantum gravity is that it must be scale invariant as well 



as reparametrization invariant [p!7| ; local shifts gab e.'^Qab of the metric scale factor are 
simply shifts of an integration variable in the functional integral over metrics. The effect 
of such a shift is to produce the trace of the stress tensor, which we can generally expand 
in terms of scaling fields 

(T.,) = = j e-^a^O' . (2.1) 



The /3i are the beta functions of the theory. The result ( |2.1| ) must vanish apart from 
possible contact terms with local operators. Often we are interested in gravity coupled to 
some number d of scalar matter fields 

^5 = 47 / Mn^) + R'^'^W + daX^dbX\G^^[X)g'''' + B^^iXY'^y^ . (2.2) 



in which case the vanishing trace of the stress tensor enforces a set of conditions on the 
spacetime fields which are the equations of motion of string theory. One may adopt either 
of two points of view: In 'critical' string theory, one sets ^-Z^atter = 0; then one has an 
additional symmetry - Weyl invariance - which one can quotient by, so that one integrates 
over the gravitational measure j^jf x w^i ' ^^^^ regards all components of the 2d 

metric as gauge degrees of freedom. On the other hand, from the 'noncritical' string point of 
view, one can trivially achieve scale invariance by integrating over all possible scale factors 
^^Diff ^ • This is equivalent to a particular class of solutions to critical string theory if we 
regard the local scale factor as another scalar field like the X's|T^; perhaps it is completely 



equivalent if we complexify the space of metrics. In either case the object is to compute 
and solve the equations of vanishing stress tensor. There are two common methods of 
calculation: the spacetime weak field expansion and the 2d loop expansion [^U| . In 



the former one expands around a known fixed point {e.g. free field theory in d spacetime 
dimensions) 

S=i^ j dXdX + 5T{X) + dX^BX^iG^JDiX) + SG^^iX) + 5B^,{X)) . 

Bringing down powers of the composite operator perturbations will reliably find nearby 
fixed points (unless one has chosen a singular parametrization of the coupling space) . The 
origin of the beta function is the overall scale divergence 

The similarity of this expression to the Koba-Nielsen prescription for the calculation of 
the string S-matrix is not coincidental [|5T|. Indeed, the effective action of string theory 
(actually any field theory) is determined from the S-matrix by subtracting the contributions 
of intermediate on-shell poles. In the Koba-Nielsen formula these are the logarithmic 
subdivergences in the integrations over the locations of composite operators; the beta 
functions are the equations of motion following from this effective action. Note that the 
kinetic operator d^S/dgidgj = d(3i/dgj for small fiuctuations is the anomalous dimension 
operator - the linearized scaling operator Lq + Lq — 2. The disadvantage of the spacetime 
weak field expansion is that one must work in a particular coordinate basis in 2d field 
space, so spacetime general coordinate invariance (invariance under 2d field redefinitions 
of the X's) is not manifest. This is an advantage of the 2d loop expansion; rather than 
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working in a specific basis of tensor fields on tlie spacetime manifold, one expands in small 



fluctuations of the coordinates X, which can be made manifestly generally covariant ||2Cl|| . 
This perturbation series is an expansion in spacetime variation of the string coordinates 
relative to the 2d Planck constant, the inverse string tension a' . A superposition of these 
two methods yields the beta functions [pO| 

=R^^ - 2V^V,I) + V^TV^T = 
(3^ = — ^ + + AiVDf - 4V^D + [VTf + ViT) = (2.3) 
P'^ =- 2\/^T + AVDVT + V'{T) = 

in a double expansion in fleld strength and a'; V{T) is a generic potential V{T) = ^T^ + 
0{T^). Both expansions are required here; the loop expansion for general covariance, and 
the weak fleld expansion to incorporate the tachyon. 

From the critical string viewpoint, we are interested in strings in d spacetime dimen- 
sions. The equations ( p.3|) are solved by {(p = X'^, say) 



D (2.4) 



Although each equation of motion in (|2.3|) is satisfled at its leading nontrivial order in 



powers of T, none is solved at subleading order. However we are only considering the 
lowest order equations; one might hope that higher orders correct the problem. After all, 
the a' (loop) expansion is not valid here since 7 is not small; the kinetic term in the tachyon 
beta function is only found after a resummation of loops, but then we have no reason to 
ignore terms involving the gradient of the dilaton fleld. Fortunately we are looking for a 
solution with rather special properties: it is ( |2.4| ) at lowest order and its renormalization 
involves only (p-dependent flelds. Since V<^T oc T, the exact solutions G, D, T are power 
series in T (assuming the weak fleld expansion is summable). Therefore we can flnd a fleld 
redeflnition - the reversion of the power series G{T), D{T), T{T) - such that (|2.4| ) is the 
exact solution0. The importance of VT cx T is that the fleld redeflnition required is local 
in spacetime. 



^ This argument is due to T. Banks. 
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The noncritical string viewpoint is somewhat more helpful in this situation. Here we 
have d — 1 string coordinates coupled to a dynamical metric gab = e'^^fjab where gab is a 
fixed metric (up to moduli). For d — 1 free scalar fields we have 

Integrating dSeff/Scj) = (3 gives an effective classical action, the Liouville action 

Self = ^ / + ^^R'^' + ^e-^^ (2.5) 

for the dynamics of the metric. The integral over metrics restores scale invariance if we can 
succeed in correctly performing the functional integral or otherwise quantizing the theory. 
From critical string considerations, one might expect a problem since the tachyon stress- 
energy appears to feed into the metric beta function. There is in fact a renormalization 
prescription such that S^ff is conformally invariant [ITtH . Since the metric is dynamical, the 
coordinate regularization scale and the covariant scale are related by a fiuctuating field: 
5 = e'^'^/^e. Therefore a regulated action should be, e.g. 

and we should calculate SSe/f/Scf) = (3 incorporating the effects of the regulator. In the 
loop expansion, normal ordering an exponential means resumming self-contractions; this 
introduces cutoff-dependence, and therefore (^-dependence. One finds that the entire effect 
of imposing (3ef / = is to renormalize the parameters in the Liouville lagrangian: fj, — > fXren 
(which can be absorbed in a shift of ^), and Q Qren = f + 7[11[0- Essentially we are 
allowed to "normal order" the Liouville term, and provided its quantum scale dimension 
is (1,1) and 7 < ^ all is well. Physically, at short distances {(p —00) the exponential 
potential term is small, so divergences should be those of free field theory. It seems that 
the complicated field redefinition of the critical string theory amounts to the change from 
coordinate to covariant regulator on the world sheet; recall that spacetime field redefinitions 
are equivalent to a change of renormalization scheme in the 2d sigma model. It is intriguing 
that in this low-dimensional context there can be nontrivial solutions of the string equations 
of motion whose stress-energy does not back-react on the spacetime metric. 

One indication of the conformal invariance of quantum Liouville theory is the degree 
of singularity of the operator product expansion of the Liouville perturbation; assuming 
free field operator products (but see [|TT| [pi}[) 

g7^W g7^(«^) _ |^_^|-27'e27<^(^) . (2.6) 
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In the weak coupling regime 7 -C 1 the singularity is integrable and the beta function cal- 
culated above will continue to vanish at first order in the perturbed theory. This reasoning 
does not give the correct upper bound on 7 however, indicating that problems set in for 
7 > 1 (c > —2) rather than the observed upper bound 7 = \f2 (c = 1). Ordinarily in 
conformal field theory, an operator that produces the identity in its operator product with 
itself has 

0{z)0{yS) ~ \z-w\-^H 

from which one can read off the dimension of that operator. It seems to be that, even 
though it does not produce the identity, in its self-product the singularity of e^'^ is four 
times its 'effective dimension'. Then from ( |2.6|) we would conclude that e"''^ becomes 
effectively irrelevant at c = 1, or more generally Oe°"^ becomes effectively irrelevant for 
OL > Q/2. This behavior is rather similar to what happens when attempting to incorporate 
massive string states in the sigma model lagrangian of the critical string. The on-shell 
vertex operator for such a perturbation has h = 1, but only by virtue of an irrelevant spatial 
operator coupled to a negative dimension temporal operator. The disease of irrelevant 
operators, namely nonrenormalizability, shows up in the appearance of many low dimension 
operators with highly singular coefficients in the self-OPE. Such a similar explosion occurs 
for operators Oe°"^ with a > Q/2 (for instance the 'black hole perturbation' p^] dXdXe^'^ , 
hence it is in no sense a 'perturbation' of fiat spacetime Liouville theory). 
To recapitulate, we have 

Tzz =0 (Liouville e.o.m.) 

These look like free field results, but it must be stressed that ((/)^) is not the free field 
propagator. The Liouville equation is geometrically the equation for constant (negative 
for /X > 0) curvature surfaces R = —fx. The semiclassical expansion is valid for 7 << 1; 
rescaling (j) — > ^(j) puts an overall factor of 1/7^ in front of ( p.5|) . If we plot the potential 
for the zero mode 
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we see that a stable solution exists only for surfaces of genus g > 1 {in the absence of point 
sources of curvature; see below). The stable point {(j)) increases with genus; in fact there 
is a scaling relation^: let (j) (j) + -loga, then 



Z{gstr, ^J^) = Z{gstr,a^)a ^^^^^ 
implying 

Zg{g,tr.lA=Cg ^i^^-^'^^ . (2.7) 

The full classical solution to the Liouville equation of motion is discovered through 
its geometrical role; e^'^ must be a density under coordinate transformations, e^'^^^^ = 
|^|2g7<?^'(2 ) Since the standard constant negative curvature metric on the upper half 
plane (UHP) is ds^ = dzdz / {Im z)^ , (f) must look locally like 



-log 
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16 dA{z) dA*{z) 
7 iA-A*y 



where A, A* are local coordinates on the surface; i.e. A{z) is the map from the UHP to 
the Riemann surface S 

4^e^^dzdz "^^"^^^ 



16- — {A-A*y ■ 

The line element on the UHP is invariant under SL(2,R) transformations A — > = 
g{A). This transformation must leave (p invariant, but may do so in a nontrivial way, e.g. 
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by making a circuit of a nontrivial closed path on S (see fig. 2). That is, the monodromy 
of the local coordinate A is a set of SL(2,1R) transformations which cover the surface S 
onto the UHP. There are several classes of monodromies: 



(1) elliptic monodromy - g = 




is conjugate to a rotation, i.e. there exists h such 



that hgh~^ = y^^^Q ~cos"/)" '^^^ surface S has a conical singularity of deficit angle 
9. Deficit angles 9 = j E 'E, are 'nice' since they are covered by the UHP; j = ^ 
requires an n-fold branched cover of the UHP at the fixed point of the rotation. 

(2) parabolic monodromy - h exists such that hgh~^ = (g i)) a translation (in terms of 
deficit angles, 9 — n and the surface has a cusp). 

(3) hyperbolic monodromy - hgh~^ = (^^ ^^^^ is a dilation. The identification of the 
UHP under g makes a handle. 

A complete solution of the Liouville equation consists of representing the surface S 
by its fundamental group, the discrete subgroup (Fuchsian group) of SL(2,1R) that covers 
the surface onto the UHP; then one constructs the automorphic function A{z) convariant 
under this group. 

The three types of monodromy above are continuously related, as can be seen by 
pinching a handle on S: 
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In the classical theory, a deficit angle ^ is a delta-function source of curvature: 



Integrating both sides, one sees that there is a solution whenever 2 — 2g + Oi/'j < 0. Eq. 
is the saddle point of the functional integral 

V(l)e-^^^^Y[e^^^''K (2.9) 

i 

On the other hand, there is no local source for hyperbolic monodromy (the fixed points of 
the SL(2, IR) transformation g on the UHP are not on the surface E). In terms of Liouville 
dynamics, the initial field configuration for (p (along some closed loop generating the mon- 
odromy) never propagates to ^ = — oo where the field configuration can be interpreted as 
localized at a point in the covariant metric. Note also that there is no classical geometrical 
interpretation for deficit angle 9 > n; i.e. when the source contributes half the curvature 
contribution -^Q j R oi the sphere in the functional integral. Two parabolic cusps turn 
a sphere into a cylinder; we cannot go beyond this while maintaining the geometrical in- 
terpretation of (j). Quantum mechanically this means |T^ — %~'^~'^'2' This 



does not mean that operators with such Liouville charge 'don't exist', rather merely that 
we cannot give them a geometrical interpretation. From the discussion of section (2), we 
conclude that gravitationally dressed operators with a> Q/2 are 'effectively irrelevant'. If 
we perturb the action by them, new dimension one operators will have to be added to the 
action to subtract singularities, a procedure that rapidly snowballs; as vertex operators 
they cannot be renormalized simply by normal ordering because the loop expansion is not 
well-behaved. It is doubtful that KPZ scaling (see below) can be maintained. 

The above considerations motivate a brief review of SL(2,IR) representation the- 
ory |2^]. Representations are labelled by their values of the quadratic Casimir C2 = j{j — I) 
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and = m + Eq, m G ^. All representations can formally be built from the two- 
dimensional representation (Zl) ^ (gfc)) = {"cd) {Zl)- Realizing the SL(2,1R) al- 
gebra on differential operators = -^Wid2, J~ = -^W2di, = ^{10182 — W2di), the 
monomial Nj^w^w^^ = Nm{wiW2)-' {wi/w2)^°^^ transforms as part of the spin j represen- 
tation, where is a normalization. There are several cases: 

(1) The trivial representation j = Eq = m = 0. 

(2) The finite dimensional representations —2j G IN, Eq = 0, m = — j, j; these are not 
unitary. 

(3) Discrete series representations V±: For "D^, we have —J + Eq = 0, —2j ^ IN, J^ + Eq = 
0, -1, -2, for V_, j + Eo = 0, -2j ^ IN, J3 + = 0, 1, 2, .... These are unitary if 
E'o G IR and j > 0. They are related to the elliptic monodromy conjugacy classes of 
SL(2,1R) and have C2 > 0. 

(4) Continuous series representations: — | < Re i?o < |, —J+Eq 7^ 0, ±1, ±2, J3+E0 = 
0, ±1,±2, .... These are unitary if i?o G M, G ilR; orEoJ G IR, < l-\Eo\. 
These are related to hyperbolic conjugacy classes in SL(2, IR) and have C2 < 0. 
How are these representations related to Liouville theory? Recall that classically 

T,^ = + From this follows[||2j 



{d^ + ^T{z))e-^^ = , (2.10) 



and from the classical solution we have 

{A -A' 




m=l,2 




(2.11) 



for j = — ^, where 



- / 1 - / ^ 

wi = ip_i 1 = —= ; W2 = W-i ,-1 



are the two solutions to ( 2.10 ). The expressions ( 2.11 ) can be regarded as a classical 
version of conformal field theory where each measurement is a sum of holomorphic times 
antiholomorphic 'chiral vertices', glued together in a monodromy invariant way to make 
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physical single- valued expressions. From the two basic solutions one can in principal build 
all exponentials corresponding to finite-dimensional representations, e.g. 

V 

m=-3,...,3 

An important issue is whether such expressions generalize to other classes of representations 
for which the sum over m is infinite. In principle one should be able to extract the answer 
from the classical solution A[z) to ( p.8|) , the saddle point of the correlation (|2.9|) ; this is 
in turn determined by its monodromy F G SL(2,1R). 

The quantum theory may have a structure similar to that of known conformal quantum 
field theories: the correlation functions would be (generically infinite) sums of holomorphic 
times antiholomorphic conformal blocks glued together to make a monodromy invariant 
object [p5[ 

a 

For instance, in current algebra conformal field theories the conformal blocks JF^, carry 
two sets of indices, an 'external' index (like i in ( p.l2|) ) which is a representation label 
for the current algebra, and an 'internal' index which is a quantum group index, since 



the monodromy acts on the blocks like the i?-matrix of a quantum group [|2^. In the 
case of Liouville, there is no 'external' symmetry group other than Virasoro, and the 
index is continuous; however there is an 'internal' symmetry group SLq{2) which is, in 
some sense yet to be understood, the 'quantization' of the classical monodromy action on 
A(2;)[0|]2^. In this respect Liouville is like other coset conformal field theories (Liouville 
can be thought of as the coset conformal field theory based on SL{2,R)/N, where is 
the Borel subgroup ||2^ ) , where the external symmetry is gauged and the internal quantum 
group symmetry remains but is 'confined' - only appearing when one pulls the theory apart 
into holomorphic constituents. In terms of the parameter 7 in the Liouville lagrangian, the 
quantum group parameter is q = exp[z7r7^] PB| p. The translation into SL(2,1R) notation 

of the condition a < Q/2 means that e^^'^ must have j < ^ = ^ + ^; the quantum 

2 77 

dimension is hj = —j — \ — !)• In these expressions the first term is the classical value; 
the second is the quantum correction. Note that Lq is also real for j = ^ -|- iA, A G IR, 
for which hj = ^[(^)2 + A2]. 

An analogy with the SU(2) WZW theory might be helpful here. The classical solution 
to the theory is 

9ab{z,z) = Ygaciz)gcb{z) 
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where the holomorphic part g{z) transforms under left multiphcation by the loop group 
of SU(2) and the Virasoro algebra, and also by right multiplication under SU(2). The 
holomorphic constituents are determined by their monodromy g(z) — * g{z) ■ h around 
nontrivial closed loops on E. gab{z) intertwines with gcd{z') via a classical r- matrix; upon 
quantization this internal index labels an SUq{2) quantum group representation with q = 
exp(-^^) ||26|| , and the intertwining is via a quantum i?-matrix. For SL(2) and for Liouville 
a major complication is that the fusion of representations closes on the continuous series 
representations and it is not clear that one has a useful or effective way of decomposing 
the correlations on intermediate states. 

Although the above analysis is rather appealing as it places the Liouville theory within 
the customary conceptual framework of conformal field theory, it has not proven to be 
sufficiently powerful to enable the calculation of correlation functions. It is a useful route 
in the case of rational conformal field theories because the monodromy representations 
are finite, related to a closed system of differential equations that one may derive from 
various Ward identities [^5| . In the case of Liouville theory, and also noncompact current 
algebra, the monodromy representations are infinite-dimensional; and the Ward identities 
are not sufficiently powerful to give a closed system of differential equations that will 
determine the correlation functions. Although finite dimensional representations e~^'^'^^^^ 
lead to finite order differential equations in z, general matter operators require infinite 
dimensional representations which don't satisfy any simple equations; such representations 
are required in intermediate states and therefore upon factorization, hence the dependence 
of correlations on the position of these operators is an open question. It seems that at least 
in the case of Liouville theory, other analytic techniques are available (see the lectures of 
D. Kutasov at this school, and references therein). 

The above analysis points to a picture of 'quantum Riemann surfaces' where the 
classical solution A{z) of Liouville theory is deformed into some kind of quantum conformal 
block, its classical SL(2, IR) monodromy deforming into an SLq(2) quantum group structure 
with q related to the coupling constant 7 of Liouville theory. 



3. Scaly behavior 

Regardless of our understanding of how to compute Liouville correlation functions, 
there are several general statements we can make about those correlations assuming the 
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quantization preserves free-field ultraviolet behavior even in the presence of the exponen- 
tial interaction term. The strongest 'theoretical' reason for such a presumption (made 
implicitly by the authors of is that physical short distance as measured in the metric 
ds"^ = e^^dz^ is where the Liouville potential is exponentially small, and so should not 
affect the free field results. The 'experimental' evidence for the validity of this assumption 
is the agreement of predicted scaling relations!]^, and more recently tree- level S- matrix 
elements, with the continuum limit of discretized random surfaces in the matrix mo del j^. 

What are the Liouville predictions? Suppose we wish to study conformal matter with 
central charge c(= d above) coupled to 2d gravity. The matter theory will contain a set of 
scaling operators (9^"**'='' of scale dimension hk- A generally covariant theory must make 
Qmatter jj^^Q ^ coordluate density of scale dimension = (1, 1); for instance, 

Qgrav _ ^CUf^ (f> ^matter -j^^ 

with 

K.tai = -^Oikiak - Q) + hk = I ■ (3.2) 

Then the integrated correlation functions of O^''"" are 2d general coordinate invariant. The 
exponential Liouville dependence of such operators will modify KPZ scaling in correlation 
functions, cf ( |2.9| ). The string theory interpretation of this gravitational dressing is that the 
Liouville field (p is the (Euclidean) time component of the string's position in spacetime; 
the time ((^) dependence of (|3.1|) is just that of a linearized solution of the (Euclidean) 
spacetime equations of motion with energy ak = jkl'i and the relation — | (a^ — Q)^ + hk = 
between ak and hk is the mass shell condition for linear perturbations, i.e. hk is the 
eigenvalue of the spatial Laplacian Lq'"""'^. This makes it clear that conformal matter 
is from the spacetime point of view the special class of stationary solutions of the string 
equations. From the previous analysis of Liouville theory, if 1 — /i^ > ^ then there does 
not exist a local generally covariant measurement corresponding to the matter operator 
Ok] i.e. O'^"-"""^ looks local in coordinates {z,z), but its gravitational dressing requires 
hyperbolic monodromy jk = + i^ which has no local interpretation in the covariant 
theory. Seibergfl^ has dubbed such operators tachyonic since their spacetime mass shell 
condition (|3.2|) implies imaginary mass for the corresponding string state in spacetime. 
The reason that some perfectly sensible local operators may not have a local gravitational 
dressing is that each measurement perturbs the local geometry by making a small deficit 
angle ak in the surface at the point of measurement; there are no probes 'outside' the 
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universe that can make such a measurement without perturbing the geometry. It can and 
does happen that when gravity is switched on the geometry is perturbed too violently by 
some operators to have a good, local continuum limit. 

In conclusion, we can take away the following main lessons about Liouville theory 
coupled to conformal matter: 

(1) KPZ scaling: 

(2) Operator scaling dimensions O^""" = e"'='^Cj^""'^'~ that shift the KPZ scaling rela- 
tion by ak in correlation functions, —^{ak — + hk = 

(3) A 'phase transition' at d = 1 where the gravitationally dressed identity operator 
becomes tachyonic in the sense described above. There is no obvious reason why tree 
amplitudes might not be analytically continued as in (i = 26 to yield a sensible classical 
string S-matrix, but (as ai d = 26) loop amplitudes are infinite because tachyons cause 
infrared divergences in loop integrals even in Euclidean spacetime. As emphasized by 
Seiberg[^, this phase transition is not always at c = 1, but occurs whenever the spectrum 
has physical tachyons hmin < 

Properties (1) and (2) are special to gravitationally dressed conformal matter; (3) 
is expected to be a generic feature persisting even when the matter theory is massive. 
One advantage of the matrix model, to which we turn next, is the ability to calculate the 
partition function and correlations even for massive matter. 

The KPZ scaling relations suggest that a large part of Liouville dynamics is accounted 
for by the zero modes. Also c < 1 minimal models coupled to gravity have the tachyon as 
the only physical state (there are some additional physical states at nonstandard values 
of the ghost number ||29|| whose meaning is less clear); reparametrization invariance cancels 
the fluctuations of the longitudinal modes, leaving only the center of mass motion of the 
string when the spacetime is 1+1 dimensional. Indeed, it has been shown that free string 
propagation is accurately described by the quantum mechanics of the zero mo des 0i 
One can imagine solving the reparametrization invariance constraints Tqq = Tqi =0 (or 
equivalently the BRST invariance constraints in conformal gauge) to show that physical 
states contain no excitations of the string's nonzero modes; the remaining constraint is 



^ String interactions are not saturated by the zero modes, which has been interpreted in [12| 
as due to the violation of the single string physical state conditions by contact interactions in the 
vertices. 
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the Wheeler-deWitt equation Tqq''° """"^^ = on the zero modes. In a conformal matter 
theory coupled to gravity this separates into the dynamics of the matter zero modes, whose 
spectrum of scahng dimensions couples to the Liouville zero mode equation pi]| [PT| 



(i—f + A^if + u'' ^'o(£) = 



with i = e^'^l'^ and v 




y). The appropriate solutions to this equation are the 



modified Bessel functions 



(3.3) 



Physically, the function K^^ = sin z^ti- i-^-'^ " ^>^] linear combination of 'incoming' 

and 'outgoing' waves I±^{2y/JIi) which is exponentially damped in the infrared £ — oo, 
indicating total reflection. Some of the strongest evidence for the viability of the conformal 
field theory approach to 2d gravity is the appearance of these wavefunctions within the 
matrix model. 

4. The matrix model 

The Hilbert space of Liouville is functionally infinite dimensional, however the sub- 
space of physical states is much smaller - at most a countable number for c < 1. Naively 
each gauge invariance removes one canonical pair of variables, one by a choice of gauge 
and another by the gauge constraint {i.e. that the generator of gauge transformations an- 
nihilate the physical subspace - the analogue of Gauss' law in QED). In 2d gravity there 
are two reparametrization degrees of freedom ^"(0- The time components of the 

metric are the Lagrange multipliers of the gauge constraints Tqq = Tqi = 0; the canonical 
degrees of freedom are the spatial metric e^'^ and its conjugate momentum tt^. Thus we 
have the possibility to gauge away up to one scalar matter fields' worth of local degrees 
of freedomi. As usual in quantum gravity the difficulty in imposing the constraints lies in 
the Hamiltonian constraint Tqq = 0, since this involves the way that the 2d worldsheet is 
carved up into spacelike hypersurfaces and therefore involves the (coordinate) time devel- 
opment. This makes it difficult to find nice global gauge invariant states. Nevertheless, the 

^ Note that in this counting the Ising model has d = ^ because tt^ = tp is 'half a canonical 



pair. 
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lesson to be drawn is that it should pay to reformulate the path integral on the space of 
physical configurations because it is expected to be much smaller than the space of metrics 
plus matter field configurations. The simplest way to enumerate physical configuarations 
is to discretize the 2d surfaces. Consider a 2d surface built by gluing together uniform 
squares each with sides of length e (or triangles, pentagons, etc. - it turns out not to 
matter which [1^ unless one artificially tunes couplings PO]). 



fig. 4 Patch of discretized surface. 

This lattice spacing replaces the covariant cutoff of Liouville theory. The only local 
freedom in pure gravity resides in how many squares meet at each vertex, which determines 
the local curvature discretized in units of 7r/2 (see fig. 4). Note that each surface is dual to 
a graph (see fig. 4) , each square being dual to a $^ vertex, each side of the square being 
dual to a 'propagator' of the $^ Feynman graph. Thus counting all graphs with A vertices 
counts all surfaces with area Ae'^; we can call this the partition function for discrete 2d 
Euclidean quantum gravity. The continuum limit consists of taking e ^ 0, A ^ oo such 
that the physical area Ap,^y^ = Ae^ is finite, assuming such a limit exists. In other words 
we concentrate on surfaces with a very large number of triangles; the statistics of these 
surfaces is governed by the large order asymptotics of graphical perturbation theory. We 
expect to be able to generate any local curvature in the continuum by coarse-graining over 
a large number of 3-, 4-, and 5+-coordinated vertices on the dual tesselation, of positive, 
zero and negative curvature, respectively. 

In the pure gravity case, the generating function for the $^ graphs dual to the dis- 
cretization is the integral 

I ci$e-^*'+»*' , (4.1) 
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i.e. the coefficient of in the asymptotic expansion at small g is the number of surfaces 
with area A in lattice units, since this asymptotic expansion is the enumeration of $^ 
Feynman graphs. The coupling g is to be identified with the bare 2d cosmological constant 
g = exp[— /Xb^^e]- Note that g is positive in order that each surface is counted with positive 
weight, so the generating function (]4.1D diverges badly; we can interpret this through 



the asymptotic expansion as the statement that the entropy of large surfaces diverges 
uncontrollably. To cut down on this entropy we can try to count only surfaces of fixed 
genus, hoping that although the sum over topology is infinite, the individual terms in the 
series might be finite. This turns out to be the case for d < 1 matter. 



Matter can be incorporated by introducing a label set for the dummy variable $||31 
then configurations are enumerated not only by the connectivity of the graph but also the 
element of the label set (which we think of as the value(s) of scalar field(s) or of discrete 
spin variables) at each point on the graph. For instance in the Ising model one considers 
two matrices M, with integrand 

exp(-tr[a6M2 + ^N^ + cMN + gdM^ + §A^^]) . 

The couplings a, b can be set to unity by a rescaling o( M, N (however such redundant 
couplings can have physical effects through contact terms in loop correlations [p2| ) . That 
leaves c, related to the Ising temperature because it controls the probability of transitions 
between up-spin (M) and down-spin (A^) subgraph domains in the diagrammatic expan- 
sion; d is related to the magnetic field since it preferentially weights one of the two spins; 
and g is our friend the surface cosmological constant. The trick in this and similar cases is 
to take the continuum limit cleverly so that both graph connectivity and label fluctuations 
approach criticality. Fine-tuning more complicated potentials yields a tower of critical 
points |]33[]. 

To count the surfaces of fixed genus consider the NxN Hermitian matrix field 
and generating function 

g-tr[i$2-(g/Ar)$4] _ 

The reason to make $ Hermitian is that the Feynman graphs carry an orientation and so 
the partition function counts orientable surfaces only. For the generalization to unoriented 



surfaces, see [^. Looking at simple graphs 
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shows us that those which can be laid out smoothly on a surface of genus g come in the 
generating function with a factor A^^~^^. Each closed index loop traces over the indices in 
that loop and gives a factor A^; the coupling constant of the generating function is chosen 
to he g/N so that adding a square without changing the topology keeps the counting 
fixed (the extra cancels the additional index trace). We have 



(4.2) 



We have the right to call this generating function the partition function of (discretized) 
string field theory in this low- dimensional situation since it is indeed the object whose 
asymptotic expansion in is the sum over surfaces of some 2d field theory describing the 
string background. The limit N ^ oo picks out the sphere contribution Cq, the classical 
limit of the associated string theory; we can evaluate this from the integral ( [4.2|) by saddle 
point techniques even though the full integral doesn't make sense (the saddle is not a global 
minimum) since the leading contribution is just the value of the integrand at the saddle. 
To find the saddle, decompose the matrix in terms of 'matrix polar coordinates' as 

$ = UAU-^ 



where U is a unitary matrix and A is the diagonal matrix of eigenvalues of $. The partition 
function becomes 



Z 



sft 



1 U-^dU j dA 






exp 


d{U,A) 





N 



1=1 



The Jacobian is easily evaluated by noting that a) it vanishes whenever any two eigenvalues 
coincide; b) it is symmetric under permutations of the eigenvalues; and c) it must scale like 
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unique function with these properties is ni<j(-^i ~ -^i)^- The justification 
for a) is analogous to the vanishing of the Jacobian from Cartesian to polar coordinates 
in ordinary multidimensional integrals: the origin is a fixed point of the symmetry group 
of rotations. The transformation is singular there so the measure must vanish. Similarly, 
the coincidence of two eigenvalues is a fixed point for an SU(2) subgroup of U(A^), so the 
Jacobian is singular. The partition function becomes 



^sft 



HdX 



-Af^^\/(A,/V]V)+2^^^.1og|(A,-A,)/y]V| 



The mental picture to adopt is to consider the eigenvalues as a set of particles lying in 
a metastable well, interacting through a logarithmically repulsive 'Coulomb' force. As 
N oo we can replace the set {A^} by an eigenvalue density (which we might regard as 
a collective coordinate or string field 0) and the saddle point equation is 

J -a IJ- J -a 



Consideration of analytic properties provides the solution 



Near the endpoint a of the distribution, p(A) behaves as p(A) ~ Va^ — A^; but when 
9 ^ 9c = the analytic behavior changes to p(A) ~ (a^ — A^)"^/^. What is happening? 
The edge of the eigenvalue density becomes softer because the last eigenvalue approaches 
an instability where the Coulomb repulsion of the rest of the eigenvalues overcomes the 
external potential force and pushes it out of the metastable well. This instability refiects a 
phase transition in the surface dynamics on the sphere: recall that the entropy of surface 
configurations grows with the number of plaquettes, but there is an energy cost exp[— /xA]. 
At g ^ Qc = e~^= from above these balance, and Z is dominated by surfaces controUably 
large compared to the cutoff e. Inserting p(A) into the saddle point action gives 



]\J2 '^saddle 



\{X)V{X)- 1 1 p{^)p{X)\og\X-^\ 



:^(a2-l)(9-a2)-iloga2 

x5/2 



'Co{9-9cj 

19 



where e^fXren = 9 — dc- This saddle point action is the partition function of 2d gravity on 
the sphere, S.^^ie = Zl'^h.r. ~ CqN'^[(? iXrenfl'^ . Letting 



9str 



gives Z^^y^^^^ = Co{iJ,rin/gstr) = Cqh^/"^] in general the full string partition function de- 
pends on the couplings iiren and Qstr only through the ratio /U. Note that to get a sensible 
continuum result we must let scale with the surface cutoff 5 ~ e^'^^'^, i.e. Qstr is dy- 
namical - just what KPZ scaling predicts! Plugging numbers into the Liouville formulae 
0)5(13) fi^d Q = 5v^, 7 = 2/v^ for (i = 0, hence Liouville theory predicts 

as observed! 

Note that the important feature of the phase transition used to take the continuum 
limit is the quadratic maximum in the effective eigenvalue potential (linear vanishing of the 
force on the last eigenvalue at (7 = g^, so it shouldn't matter whether we use $^ squares 
or triangles in the microscopic theory. In other words, only the rate of vanishing of p 
near its endpoint is universal; the details of the eigenvalue distribution far from this region 
as well as the value of gc are irrelevant quantities. By fine tuning the relative weight of 



different polygonal simplices one can however reach other sorts of critical points ||30|] [p6 |. 

The advance embodied in the discrete approach is the ease with which one may cal- 
culate the coefficients Cg in the partition function[0, as well as all correlation functions. 
To date it has only been possible in the Liouville approach to calculate (after much effort) 
correlation functions on the sphere. To go beyond the sphere one needs to among other 
things incorporate corrections to p(A) due to the discreteness of eigenvalues; p is not a 
smooth function at the level but rather a sum of delta functions. In fact the whole 
methodology above is rather cumbersome, and it proves simpler to reformulate the prob- 
lem. It is apparent that the main difficulty is the Coulomb interaction among eigenvalues, 
so we might try to find variables that diagonalize this interaction as much as possible. This 
happens when one writes the square root of the Jacobian as a Slater determinant 

n(A.-A,)=det,,[V'.(A,)] = |vl>(A)) , 
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where ipi{X) = A* — 1 + Ci_2A*~^ + . . . + cq. The q are arbitrary since they correspond 
to — > V'i + V'fe) k < i, which cancels from the determinant. A convenient choice is to 
orthogonahze ip with respect to the measure e" 



*t(A)) = (^t(A) *t(A) 



Letting 



we find 



hnSnm = J V'n(A)e ^^^Vm(A) , 



(4.3) 



N 
i=l 



up to a t-indcpcndcnt normalization. The computation of the partition function reduces 
to the determination of the hi{i). To this end, define 



h„Pnm = J V'n(A)( 

Considering different matrix elements yields 



-V{X) 



d 

dX 



V'm(A) 



(4.4) 
(4.5) 



\ 



1 62 a2 

1 63 as 

Vo ••• ••• 



(4.6) 



This also tells us (by acting A to the left) that = hi+i/hi; a little thought shows 6^ = if 
V is an even function (since then we have the symmetry V'j(A) = (— )-'V'(~A)). Integration 
by parts for P gives 

P = V'{Q)+ , (4.7) 

where the + subscript means taking only the upper triangular part, setting the lower 
triangular and diagonal entries to zero, since d/dX must lower the order of similarly a 
— subscript refers to the lower triangular part. Now 



[^,A] = 1^[P,Q] = 1 
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(4.8) 



is an equation for a, b which solves our problem, i.e. 

N N 



Zsft = Whi = hQ exp[^(A^ - i)\ogai 



i=l i=l 

AfRcionados of integrable systems will recognize that iiV = ^A^ then P, Q is the Lax pair 
of the Toda chain; actually trjQ'^"'"^} is the /c*^ conserved charge of the Toda problem and 
(Q'^)+ is the k^^ Toda Hamiltonian. In other words, 

^ = [Qi,Q] (4.9) 



are commuting flows. One easily checks that these flows preserve the string equation ( ^.(jj ) 
Finally one can show that Zgjt is a particular type of 'tau-function' of the Toda system; this 
conveys little information but sounds impressive, the content being primarily in equations 
( [4 .81) and ( [4.9| ). Commuting flows in this problem should not surprise us; we must have 
a? = iZ t, , etc., so of course all the derivatives commute. The nontrivial fact is that 

atidtj dtjdti ' ' 

we can represent Z as an inner product in a Hilbert space 

Zsft{t) = (^^|e-S(*^-*'=)^' (4.10) 

where the derivatives with respect to the couplings act as operators. The different flows 
in coupling space then form an integrable system in this Hilbert space. 

Again we need to take the continuum limit. To this end it is convenient to absorb 
a factor e~^/^ in normalize i(jk -^^"^k and take V even (although the latter is not 
essential). Then the Lax pair P, Q are transformed into 

Q = Q+ - Q_ P = V'iQ)+ - V'iQ)- . 

The scaling limit consists of taking const. + e'^u{x), where = 1 — e^x. The second 
order difference operator Q becomes (in the sense of matrix elements) a second order 
differential operator + u(x); ^ Zn+iZn-i yjg}(^g ^(^^.^ _ Q'ijr where Zgft = e~-^. 

Finally P scales to some odd order differential operator depending on the details of the 
potential (we made the rescaling in the wavefunctions in order to make this manifest). 
Generically ^ which is boring, i.e. 



[P, q,] = \ = u' = T'" ^T = x^ 
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which corresponds to p ~ y/a"^ — X^. Fine tuning as before should yield p ~ (a^ — A^)"^/^ 
corresponding to P — 9'^ + vi{x)d + vq. Then the equation [P, Q] = 1 is the continuum 
statement 

u + uu = 1 , 

the Painleve I equation. This equation embodies through its asymptotic expansion in 
X = /i the entire perturbation series for Zgft about this background. The leading solution 
at large x is = x, which gives ~ x^/^ which is the famous KPZ scaling relation! 
Plugging this solution into the differential equation and iterating yields 



2.1(2-23) 



3=0 



SO we see (5/7 = 5/2. Note that derivative terms in the string equation come with factors of 
the string loop coupling Qstr ~ 1/-^, which is why the genus expansion expands in powers 
of derivatives of u{x). Douglas[jl^ has generalized this construction for an arbitrary 



pair of differential operators P, Q of integer order p, q satisfying (|4.8[ ), which produces BPZ 
minimal matter coupled to 2d gravity. Changing the bare matrix potential scales to 
continuum perturbations of P, Q by lower or higher order differential operators for relevant 
and irrelevant perturbations, respectively [^@. In the two (and more) matrix case, we 
can again eliminate the angle variables in the partition function[B3[. To 'diagonalize' the 



eigenvalue interaction, choose independently the left and right wavefunctions (\l/t| and 
l^t)!^^]- Defining matrices M, Pm and A^, Pn analogous to ( |4.4D and ( |4.5D , we have the 
variational equations 

PM=N + VitiM) 
Pn=M + VI^{N) 

The number of nonzero diagonals of M, are determined by the degrees of Vn, Vm 
respectively. Tuning these potentials, we can make M, Pm scale to a pair of differential 
operators 

Q =d'' + Uq-2d''-^ + . . . + uo 

P =dP + Vp-idP-^ + . . . + Vo 

of orders {p,q)- The Heisenberg relation [Pm,M] = 1 again determines the coefficients 
Un(x,t), Vm{x,t) and hence the free energy d^J-' = Wg_2[|l^[^. These critical points 



^ Of course any fixed lattice operator is a sum over continuum operators with cutoff dependent 
coefficients. 
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have been identified with (p, q) conformal matter coupled to 2d gravity. The flows in the 
couphngs t are governed by the 2d Toda hierarchy on the lattice[0, which scales to the q- 
reduced KP hierarchy in the continuum [|16|| [P7| . The coefficients Cg ~ (2/i)! for large g, 
reflecting the divergence of the sum of the perturbation series[|l^[0. It is very interesting 
that these coefficients can be associated with an auxilliary, higher-action saddle of the 
matrix integral corresponding to the eigenvalue configuration where the last eigenvalue 
is moved to its unstable equilibrium 



which might be interpreted as a 'string instanton' mediating some kind of vacuum decay, 
although the precise meaning of this configuration remains unclear. 

An interesting and calculable set of correlation functions in these systems are corre- 
lations of the loop operator W{£) which cuts a hole of boundary length i in the surface. 
In the matrix model this operator is ytr$^, which in graphs inserts a source of I external 
lines graphically dual to a loop or hole of boundary length le. For instance, fig. 4 is a 
contribution to / = 24. In the continuum limit Q — > const. -|- e(— 9^ + u), I i/e this 
becomes the heat operator e~^^~^ "'"'^^ 

tr$' ~ (const. + e{-d^ + u)fl' ~ e-^(-^'+") = W{IL) 

up to an unimportant (nonuniversal) normalization. The calculation of loop correlation 
functions thus reduces to a set of convolutions of heat kernels [P3|, e.g. 

/oo p-^l 
dx dy{x\e-^^^\y){y\e-^^^\x). (4.11) 
-/i J — oo 

For the general (p, q) model we simply use the operator Q determined by the string equa- 
tions ( [4.8| ) . The loop correlations serve as a kind of generating function of local operator 
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correlations, in the sense that the £ — limit of W{i) can be recast as a sum (integral) 
over the spectrum of local operators [0. Geometrically, the zero mode i = § e''"^/^ is dual 
(in the sense of Fourier transform) to its conjugate Liouville momentum p = § ^4, whose 
eigenvalue is the deficit angle; hence it should be no surprise that a Green's function at 
fixed small i should be an integral transform over local curvatures. The small i expansion 
of (|4.11| ) is the asymptotic expansion of the heat kernel, which is evaluated in terms of 
formal fractional powers of (5[^[|15|[0. For instance, (|4.11|) expands as 



CO 



{W{£)W{£')) ~ Yl ^ (C'n ^^W> 

n=0 

(O^'^Wii)) = J'^ dx [Ql^^ {x\ e-^« \x)] 



For the definition of fractional powers of differential operators and their uses, see e.g. 
lUllM^lEll- The basis of operators defined by isolating individual terms on the RHS do 



not coincide with the dressed operators Oa of conformal field theory ||T^ . Instead, one 



must take linear combinations of O^^ with coefficients defined by the 'incoming wave' 
In/q{2y/lli'). Then one finds the two-point function on the sphere of a CFT scaling operator 
and a loop operator is 

(O^Wii)) = {2^rK^{2^e) . (4.12) 

Thus one can justifiably claim that the matrix and Liouville approaches coincide where cal- 
culations can be done in both. Moreover, according to the ideas of Hartle and Hawking[Q, 
one expects that the correlation function of a scaling operator Ok with a loop operator 
W{i) is the corresponding wavefunction of that operator - the operator creates the ap- 
propriate state in the Hilbert space of physical states, and W{i) is the probability of that 
state propagating on a Euclidean surface of average curvature — /i to a boundary of length 
i. This expectation is indeed borne out by the explicit calculations |jT^ outlined above. 

The matrix model is an effective tool for the calculation of correlations of integrated 
local scaling operators. While this is an important class of measurements it by no means 
exhausts the list of interesting questions one can try to ask in 2d gravity. These integrated 
correlations are what might be called extensive measurements, things like the area of 
the surface, or the magnetization or average energy density in the Ising model. Truly 
local measurements, where one tries to 'build a laboratory' on a patch of the surface and 
make measurements in it, have not been addressed yet because they are intrinsically more 
complicated. One must describe the laboratory in a coordinate invariant way, e.g. by 
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prescribing the geodesic distances and relative orientation of the objects in it. Even the 
simplest such 'local' measurement, the Hausdorff dimension of the surface, turns out to be 



very difficult to investigate p2|. It is of course possible that, like the scaling operators of 
large positive Liouville momentum, such measurements in two dimensions are plagued by 
strong fluctuations that render local questions meaningless. In such a situation the simple 
model does not retain the flavor of its higher dimensional cousin. Another interesting issue 
is whether we can deflne Minkowski 2d gravity independent of working in Minkowski string 
theory. Usually we induce the continuation to Minkowski signature world sheet through 
the analytic continuation of string amplitudes from Euclidean spacetime and its associated 
ie prescription. Nevertheless it would be odd if the beautiful physics of 1 + 1-dimensional 
scattering theory could not be put on a fluctuating geometry. Can we not discuss the 
S-matrix of the Ising model, or Potts or sine-gordon models, in deSitter 2d gravity? This 
also may hinge on how to deflne localized asymptotic states in a fluctuating background 
geometry. 
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